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ABSTRACT 

We present new mid-infrared and submillimetre observations for a sample of eight high 
redshift type-2 QSOs located in the Chandra Deep Field South. The sources are X-ray ab- 
sorbed with luminosities in excess of 10^^ erg s~^. Two of the targets have robust detections, 
S/N > 4, while a further three targets are marginally detected with S/N ^2.5. All sources 
are detected in multiple mid-infrared bands with the Spitzer Space Telescope. The multiwave- 
length spectral energy distributions (SEDs) of the type-2 QSOs are compared to those of two 
local ultraluminous galaxies (Arp220 and IR22491) in order to assess contributions from a 
star-forming component in various parts of the SED. We suggest that their submillimetre emis- 
sion is possibly due to a starburst while a large fraction of the mid-infrared energy is likely to 
originate in the obscured central quasar Using the mid-infrared and submm observations we 
derive infrared luminosities which are found to be in excess of L> 1O^^L0. The submillimetre 
(850 /im) to X-ray (2 keV) spectral indices (asx) span a wide range. About half of the type-2 
QSOs have values typical for a Compton-thick AGN with only 1 per cent of the nuclear emis- 
sion seen through scattering and, the remaining with values typical of submm-bright galaxies. 
Combining the available observational evidence we outline a possible scenario for the early 
stages of evolution of these sources. 
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1 INTRODUCTION 

The discovery that most local spheroids contain massive black 
holes (MBH) in their centres has renewed the interest in the phys- 
ical properties of active galactic nuclei. A number of correlations, 
such as the one between the masses of MBH and the properties 
of the host galaxies (e.g. Gebhardt et al. 2001, Merritt & Fer- 
rarese 2001) and the proportionality between the mass and the 
velocity dispersion of the stars, suggest a direct link between the 
formation/growth of the black hole and the stellar mass of the 
galaxy spheroid. 

The availability of major facilities in the submm (e.g. SCUBA 
on the James Clerk Maxwell Telescope and MAMBO on IRAM) 
led to the discovery of a large population of luminous infrared 
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galaxies at high redshifts (e.g. Hughes et al. 1998, Bertoldi et al. 
2000, Borys et al. 2003). FoUowup observations in CO molecu- 
lar gas-line emission suggest that such galaxies are in fact mas- 
sive (e.g. Greve et al. 2005, Chapman et al. 2008). The launch of 
the Spitzer Space Telescope (SST, Werner et al. 2004) has shed 
new light into luminous dusty high redshift galaxies. Sensitive mid- 
infrared observations have unveiled the presence of obscured AGN 
whose properties would have otherwise remained unnoticed (e.g. 
Lacy et al. 2004, Stern et al. 2005, Alonso-Herrero et al. 2006 here- 
after AH06, Donley et al. 2007). 

There are two obvious ways one can investigate the coeval 
formation of black holes and galactic bulges. One approach is to 
follow up on known mm/submm sources with X-rays. Somewhat 
surprisingly, however, only a small fraction of SMGs contain a vis- 
ible luminous QSO (e.g. Fabian et al. 1999, Almaini et al. 2003). 
Such a result implies that the epoch of black hole formation does 
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not coincide with the epoch of major star-forming activity. Using 
deep X-ray data Alexander et al. (2005) found that ~75% of the 
Chapman et al. (2005) SCUBA sources do contain an AGN how- 
ever, the bolometric output is dominated by intense star formation. 
This claim has more recently been confirmed by mid-infrared stud- 
ies of submm-luminous sources by Pope et al. (2008). The sec- 
ond approach is to select "obsciu-ed AGN" and follow them up 
in the submillimetre wavelengths. Such an approach has already 
been adopted by Page et al. (2001), Archibald et al. (2001) and, 
Stevens et al. (2005) with encouraging results. These authors se- 
lected targets from the ROSAT surveys (e.g. Page et al. 2001), 
XMM-Newton and Chandra (e.g. Page et al. 2003) that span the 
redshift range l<z<3, luminosity range 44.1<log Lx <44.4 and 
absorbing column densities 21<log N/j <23. 

The advantage of the present work is that we are using a sam- 
ple of "heavily" X-ray absorbed (log N/f > 23 cm~^) and very lu- 
minous (log Lx > 44 erg s"'^) AGN, the so called type-2 QSOs, lo- 
cated in the Chandra Deep Field South (CDF-S) and selected from 
Szokoly et al. (2004). The first submm detection of such a heavily 
absorbed QSO has been reported in Mainieri et al. (2005). 

In this paper we combine the submm observations of the en- 
tire sample of eight type-2 QSOs with observations obtained with 
SST We discuss the spectral energy distribution of the sources and 
compare them with local templates. We investigate a possible link 
between type-2 QSOs and SMGs, especially those known to con- 
tain a black hole. We explore the origin of the mid-infrared and 
submillimeter emission and finally discuss a possible evolutionary 
scheme for type-2 QSOs. Throughout the paper we assume a cos- 
mology Om = 0.3, Oa = 0.7 and Ho = 70 km s"^ Mpc "\ 



2 SAMPLE SELECTION AND OBSERVATIONS 
2.1 The type-2 QSO Sample 

The deep (~1 Msec) Chandra survey of the Chandra Deep Field 
South (CDF-S, Giacconi et al. 2002) resulted in the discovery of a 
unique sample of heavily X-ray absorbed (log Nif > 22 cm~^) and 
very luminous (log Lx > 44 erg s^^) AGN, the so-called type-2 
QSOs. The type-2 QSOs for this program have been selected from 
the deep (~ 1 Msec) survey of the CDF-S (Szokoly et al. 2004). 
From the 10 objects contained in the original list, 7 objects were 
observed with the JCMT within the allocated time. To this sample, 
source XID 901, was added. Although the source does not meet the 
formal Lx criterion dscribed above, it is a heavily absorbed object 
at 2 > 2. In Table 1 we list properties of the sample of type-2 QSOs 
studied here: names, coordinates, redshifts. X-ray luminosities and 
absorbing column densities. We note that X-ray luminosities are 
corrected for absorption (values taken from Tozzi et al. 2006). 



2.2 Submillimeter data 

Observations at 850 fim were carried out on the JCMT during 
2004 August and November. We used the Submillimetre Common 
User Bolometer Array (SCUBA, Holland et al. 1999) in photome- 
try mode, in which the source is placed on the central bolometer of 
the array and the secondary mirror is jiggled in a 3 x 3 pattern with 
2 arcsec intervals with a 1 sec integration at each position. The sec- 
ondary mirror was chopped 45 arcsec in azimuth at a frequency of 
7.8 Hz and nodded between the source and reference positions ev- 
ery 18s. We placed each source in the central bolometer (H7) and 



used the median of the remaining bolometers for additional sky re- 
moval. We used Uranus to calibrate the derived flux densities. Cali- 
bration uncertainties are about 10% at 850 jj,m. The pointing of the 
telescope was checked frequently while the sky opacity was mon- 
itored via regular skydips using the JCMT Water Vapour Monitor 
and the CSO (Caltech Submillimetre Observatory) Tau Meter. Ob- 
servations during August 2004 were made under very good weather 
conditions with the 225- GHz sky opacity as measured at the adja- 
cent CSO, r225, in the range 0.05-0.13 (at the airmass of the target). 

We reduced data using the standard STARLINK software col- 
lection SURF. After compensating for the nod the data were flat- 
fielded and corrected for atmospheric extinction. Each jiggle in 
turn was then corrected for residual sky noise which is correlated 
across the SCUBA field of view and often dominates the signal 
from faint sources. After sky-subtraction the data were clipped at 
the 3-a level. Since data for the same source were collected on dif- 
ferent nights we tested them for consistency with one another us- 
ing a Kolmogorov-Smirnov test, rejecting anything below the 5% 
mark. Submm fluxes are reported in Table 2. Two sources were de- 
tected at S/N ^ 4, while a further 3 are observed at S/N ^ 2.5 
level. We note that CDFS-202 is the only type-2 QSO for which the 
SCUBA observations yielded no signal at all and therefore a formal 
upper limit (corresponding to the sensitivity limit achieved by the 
observations) is quoted. 

2.3 Spitzer Observations 

The CDF-S was observed at 24 and 70 fim with the Multiband 
Imaging Photometer for Spitzer (MIPS, Rieke et al. 2004) and the 
InfraRed Array Camera (IRAC, Fazio et al. 2004) covering a to- 
tal area of 1 .5x0 .5 as part of the MIPS Guaranteed Time Ob- 
servations (GTO). Additional observations of the central region of 
CDF-S were taken as part of the Great Observatories Origins Deep 
Surveys (GOODS, for an overview see Dickinson et al. 2003). Fi- 
nally, additional 70 /xm data were taken as part of the FIDEL (P.I. 
M. Dickinson) program. 

The 24/im source extraction and photometry is described 
in detail in e.g. Papovich et al. (2007). In summary, all sources 
were treated as point-like, given the ~5 .8 FWHM angular res- 
olution. PSF fitting and flux measurements were performed using 
DAOPHOT packages. The resulting 5 a flux density limit was 80 
fj,Jy reaching a completion of 75% at this level. IRAC observations 
were carried out at 3.6, 4.5, 5.8 and 8.0 /um. Source detection, ex- 
traction and photometry was performed in a similar manner to the 
MIPS data i.e. using DAOPHOT routines. IRAC photometry was 
carried out using a PSF of l'.'8-2'.'0. The aperture fluxes in each 
band were subsequently corrected to total fluxes using known PSF 
growth curves from Fazio et al. 2004; Huang et al. 2004. The re- 
sulting 5a flux density limits were 1.73, 3.02, 10.96 and 8.3 /itJy at 
3.6, 4.5, 5.8 and 8.0 /jm, respectively. 

The Spitzer photometric data are listed in Table 2. Only one 
(XID 901) of the sources was detected at 70/um (B. Weiner, prov. 
comm.). 



3 SPECTRAL ENERGY DISTRIBUTIONS 

In Figure 1 we show rest-frame (Spectral Energy Distributions 
(SEDs) for the present sample of type-2 QSOs. To construct the 
rest-frame SEDs we used optical photometry from HST/ACS ob- 
servations (as part of the GOODS survey, Giavalisco et al. 2004), 
near-IR ISAAC imaging data (Retzlaff et al. 2009), IRAC/MIPS 
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Table 1. Optical coordinates, spectroscopic redshifts, X-ray luminosities and absorbing column densities 



nji RA(2000) Dec(2000) l.x[0.b - 2KeV]'^'^ Lx[2 - lOifeVl^.s N/j/ 10^2 

lO^^ergs-i 10*-*ergs-i (cm-^) 



202 


03 32 29.86 


-27 51 05.8 


19.8 


5.69 


150.0 


3.700 


54 


03 32 14.61 


-27 54 20.7 


3.03 


0.88 


10.67lt;f, 


2.561 


45 


03 32 25.68 


-27 43 05.7 


4.03 


1.10 


8.1913.02 


2.291 


263 


03 32 18.83 


-27 51 35.6 


8.67 


2.91 


150.0 


3.660 


27 


03 32 39.67 


-27 48 50.5 


5.23 


2.13 


28.081'^;^^* 


3.064 


112 


03 31 52.07 


-27 53 28.2 


7.52 


1.16 




2.940 


901 


03 32 35.78 


-27 49 16.82 


0.87 


0.13 


Ig 94+17.86 


2.578 


51 


03 32 17.16 


-27 52 20.7 


5.88 


1.02 


22.42l2.85 


1.099 



^ : IDs are from Giacconi et al. (2002) 
2 : Njj values from Tozzi et al. (2006) 
: X-ray luminosities corrected for absorption. 



Table 2. Mid-Infrared and submm fluxes and Itr uncertainties 



ID 


3.6 


4.5 


5.8 


8.0 


24 


70 


850 




(/iJy) 


(A'Jy) 


(MJy) 


(MJy) 


(MJy) 


(/^Jy) 


(mjy) 


202 


3.82(0.5) 


1.63:1 


4.69(0.5) 


10.47(0.5) 


76.10(0.5) 


0.0 


2.0(2.6) 


54 


4.4(0.5) 


0.02 


10.1(1.4) 


5.4I 


60(9) 


0.0 


9.74(3.65) 


45 


14.6(2.2) 


21.8(3.5) 


53.8(5.0) 


124.6(23) 


480(39) 


0.0 


2.22(1.09) 


263 


7.68(0.6) 


6.82(1.1) 


11.3(1.9) 


45(3.5) 


63(8) 


0.0 


4.82(1.11) 


27 


13.6(1.2) 


15.5(2.2) 


18.9(3.4) 


24.6(3.0) 


154(22) 


0.0 


4.81(3.45) 


112 


6.4(0.7) 


9.6(1.0) 


18.1(2.7) 


34.7(3.7) 


301(49) 


0.0 


9.87(1.76) 


901 


14.9(1.5) 


16.7(1.9) 


30.3(3.5) 


41.0(3.9) 


520(45) 


3320(1850) 


8.8(3.6) 


51 


66(5.7) 


73(6.5) 


0.0 


99(22) 


0.0 


0.0 


2.89(1.06) 



1 : denotes a (3 it) upper limit 
2: source completely undetected. 



Spitzer observations, and X-ray data (Giacconi et al. 2002). The ra- 
dio data come from the VLA survey of the CDF-S, Kellerman et 
al (2008), Miller et al. (2008). Spectroscopic redshifts for all the 
type-2 QSOs presented here are reported in Szokoly et al. (2004). 
The overall shape of the SED of the type-2 QSOs is rather different 
from the standard (type-1) QSO SED (see discussion for CDFS-263 
in Mainieri et al. 2005). From the shape of the SEDs presented in 
Figure 1 and the discussion in Mainieri et al. (2005) it is clear that 
type-2 QSOs and, in that instance CDFS-263, do not show the fea- 
tureless steep continuum expected for a typical QSO. Instead, the 
shape of their SEDs implies the presence of an additional starburst 
component. 

Here we adopt a slightly different approach and using 
templates chosen from local Ultra-luminous Infrared Galaxies 
(ULIRGs) we try to assess the dominant contribution of starburst 
and/or AON components in the various parts of the type-2 QSO 
SEDs. Our choice of templates (ULIRGs) was driven by the de- 
rived far-infrared luminosities of type-2 QSOs which are discussed 
in section 6. We stress that the aim of this exercise is not to find 
out whether there is an AGN in these type-2 QSOs (by definition 
they have a central massive accreting black hole), instead, we want 
to assess the relative contribution of the the two components (AGN 
or starburst) in various parts of the SED. In particular we focus on 
the near-infrared and submm wavelengths; it is clear that the mid- 
infrared is dominated by emission from the AGN. We chose two 
templates, the highly obscured ULIRG archetype Arp220 and the 
more typical (less obscured) starburst-dominated ULIRG IR22491 
so that we can explore the effect of obscuration in the type-2 QSOs. 

In Figure 2 we show the two templates and the SEDs of the 



type-2 QSOs presented in this study. The SEDs of the two tem- 
plates as well as those of the type-2 QSOs have been normalised in 
the submm. It is clear that in all cases (all type 2 QSOs) an extra 
AGN component is inmiediately needed to explain the mid-infrared 
part of the SED. Two of the sources, CDFS-27, and CDFS-263 are 
consistent with the SED of IR22491 (ie a typical ULIRG) with an 
extra component that comes out in the rest-frame 6-30 jim range. 
The near-infrared points follow the SED of the template implying 
that the near-infrared light comes mainly from the stellar popula- 
tion in the (massive) host galaxy. The SED of CDFS-45 rises more 
steeply than the SED of IR22491, possibly due to a stronger under- 
lying AGN and/or less obscuration. 

Sources CDFS-112, CDFS-901 and CDFS-202 are consistent 
with the Arp220 SED plus again an additional AGN component 
to account for emission in the mid-infrared. The difference with 
these sources is that they require a more heavily absorbed starburst 
component to match the near-infrared part of the SED. We cau- 
tion however, that this finding does not imply an abundance of Arp 
220-like objects at redshift of ~3, instead, we suggest that these 
sources are heavily obscured with rather unusual properties. Rigby 
et al. (2008) discuss the the lack of appropriate templates for tar- 
gets at redshifts 2-3. The remaining source CDFS-54, poses a real 
challenge when fitting the SED as it is undetected at 8 but detected 
at 24 microns. The SED of CDFS-54 falls below the SEDs of both 
templates considered here. 

Although the current sample is small, the agreement between 
the templates and the SEDs of the type-2 QSOs suggests that type- 
2 QSOs might have a significant starburst component. Star-forming 
processes could account for both the visible/near-infrared light as 
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Figure 1. Spectral Energy Distributions of type-2 QSOs. The photometric points include optical HST/ACS data (Giavalisco et al. (2004), near infrared data 
(Rettzlaff et al. 2009), Spitzer and SCUBA data presented here (shown in Table 2), X-ray detections from Giacconi et al. (2002) and radio data from KeUermann 
et al. (2008) and Miller et al. (2008). Errorbars are shown for the submm data and the MIPS 70 /im data. In the remaining cases the errorbars are smaller than 
the plotting symbols (see Table 2). 



well as the submm. The mid-infrared, however, is in all cases pow- 
ered by the central AGN. We note that the SEDs of the type-2 QSOs 
show a greater variety in shape than the SEDs of SCUBA galaxies 
(c.f . Pope et al. 2008). In the case of the type-2 QSO population this 
may reflect the various stages of evolution of the galaxies which 
may also affect the energetics (starburst vs central AGN). 

Interestingly, two of the QSOs presented in our study are also 
part of the mid-IR selected powerlaw sample of AH06. CDFS- 
901 and CDFS-112 have been classified as Narrow-Line AGN 
(NLAGN) and ULIRG-like, respectively. Since, however, the di- 
vision between the two categories is rather arbitrary we will, as 
in AH06, consider them as one class. According to AH06, galax- 
ies in this category show a steep powerlaw continuum with a less 
prominent 1.6 micron bump. Both CDFS-901 and CDFS-112 fol- 
low the Arp 220 template SED (Figure 2), except for the mid- 
infrared where the additional need for a steep power-law (AGN) 
is evident. Overall, we find a dichotomy in the shape of the type- 
2 QSO SEDs which likely reflects the varying contribution of the 
starburst component to the overall energy output of the sources. 
While starburst-related activity is responsible for the energy re- 



leased in the submillimetre, the mid-infrared is likely the place 
for the dominance of the central compact source which we discuss 
next. 



4 MID-INFRARED CONTINUUM 

At the redshift of our type-2 QSOs the IRAC bands probe the rest- 
frame near-infrared continuum while, the MIPS 24/im band probes 
the mid-infrared around 6/im which is significantly affected by the 
obscuring dust. The 24/L(m flux could originate either from emission 
related to star-forming processes (at this redshift the 6.2 fim Poly- 
Aromatic Hydrocarbon (PAH) feature enters the MIPS 24^m band) 
or it could contain indirect contributions from the AGN via heating 
of the surrounding dust which in turn emits in the mid-infrared. Due 
to the inherited complexity of the mid-infrared part of the spectrum, 
we will not attempt a decomposition of the mid-infrared part of the 
SED of type-2 QSOs but instead opt for simple colour-colour mid- 
infrared diagnostic diagrams. 

A number of authors have published results on AGN selec- 
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Figure 2. Comparison of the Type 2 QSO SEDs with those of template galaxies Arp220 (long dashed line) and IR22491(solid line). The template SEDs have 
been normalised at 850 /im. The photometry used to construct the SEDs of type 2 QSOs has been described in Figure 1. 



tion based on mid-infrared colours (e.g Lacy et al. 2004, Stern et 
al. 2005, Hatziminaoglou et al. 2005 but see Barmby et al. 2006, 
Rigby et al. 2008 for a discussion of the various criteria in detail). 
The Lacy et al. (2004) plots involve IRAC and MIPS bands only 
and are, thus, more appropriate for our purposes. Here, we exam- 
ine the mid-infrared colours of our sources and compare them to 
those of high-z submillimetre luminous galaxies (hereafter SMGs). 
bf Based on the shape of their SEDs, SMGs can be crudely divided 
in those that display a clear 1.6^m "stellar bump" and their over- 
all SED shape is similar to that of Arp220 (referred to as "cold" 
SMGs) and, those with a steep powerlaw SED similar to the SED 
of Mrk231 (referred to as warm, see e.g. Egami et al. 2004). Fig- 
ures 3a and 3b show the Lacy et al. (2004) mid-infrared colour- 
colour criteria for our type-2 QSOs and SMGs from Egami et al. 
(2004) and Ashby et al. (2006). In the IRAC plot (Fig. 3a) we find 
that CDFS-27, CDFS-45 and CDFS-263 display the reddest colours 
both in 58.0/54.5 and S5.S/S3.6 and thus are similar to those of 
"warm" SMGs. These three type 2 QSOs were best fit by the less 
obscured IR22491 template (of course as discussed in Section 3 an 
AGN component is necessary to explain the mid-infrared part of the 
spectrum). Similarly, type 2 QSOs following the Arp220 SED show 



IRAC colours with values closer to those of "cold" SMGs. In the 
IRAC/MIPS plot (Fig. 3b) the colours of type 2 QSOs vary over a 
wider range (especially the 24/5.8 ratio). AGN display colours that 
would fall within the shaded area in Fig. 3b. As expected, "warm" 
SMGs fall within the shaded region. Among the type 2 QSOs, those 
that were fit with the less-obscured ULIRG template IR2249I (plus 
AGN) fall within the shaded region expected for AGN. The remain- 
ing type 2 QSOs fall outside the expected range for AGN. 

Although the present type-2 QSO sample is small, the mid- 
infrared colour-colour plots allow us to draw some qualitative con- 
clusions. First, in both plots the less obscured type-2 QSOs (those 
fit by IR2249I) display colours akin to those of AGN-dominated 
or "warm" SMGs. In the IRAC plot, the 8.0/4.5 ratio provides a 
good measure of the AGN contribution (at z~3 the IRAC 8 nm 
band corresponds to K rest-frame). Warm SMGs and less obscured 
type-2 QSOs both display the highest 8.0/4.5 values as expected 
from hot-dust contribution to the rest-frame K-band emission. In 
the IRAC/MIPS plot (Fig. 3b) type-2 QSOs fit by IR2249I (less 
obscured) and "warm" SMGs show very similar colours. Those 
type-2 QSOs that follow the SED of Arp220 (ie the starburst com- 
ponent is more obscured) show the most extreme values in 24/8 
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and 5.8/3.6 with some type-2 QSOs falling outside the shaded area 
expected for AGN type objects. We attribute the wide range vari- 
ation in the S24/S5.8 ratio to the varying degree of obscuration, 
contribution from the central AGN and possibly differences in dust 
composition. 

Sturm et al. (2006) presented mid-infrared IRS spectroscopy 
of a small sample of primarily lower-redshift type-2 QSOs 
(0.2<z<1.3) but also included CDFS-202. They tried to estimate 
the AGN contribution to the mid-infrared (6/im) part of the spec- 
trum by performing a decomposition of their IRS spectra using a 
starburst (M82) and and AGN (linear continuum) template. Their 
results, together with the surprisingly absence of PAH features in 
the IRS spectra, indicate that the AGN contributes significantly 
(if not dominantly) in the mid-infrared. More recently, Martinez- 
Sansigre et al (2008) presented IRS spectra of a sample of z~2 
t5'pe-2 QSOs (their selection was based on radio and mid-infrared 
criteria, see Martinez-Sansigre et al. 2006) and concluded that al- 
though they are predominantly continuum dominated a large frac- 
tion shows deep silicate absorption features and PAHs. The strength 
of the PAHs is similar to that found in submillimetre luminous 
galaxies (e.g. Valiante et al. 2007). The PAHs are thought to origi- 
nate in the host galaxy and are indicative of vigorous star-forming 
activity. A similar conclusion has been reached by Zakamska et 
al. (2008) based on mid-infrared spectroscopy of optically selected 
t5'pe-2 QSOs. Combining our discussion here with the SED fits we 
also conclude that the mid-infrared part of the spectrum is domi- 
nated by emission from the AGN. 



5 FAR-INFRARED LUMINOSITffiS, STAR FORMATION 
RATES AND MASSES 

In Table 3 we present far-infrared luminosities (Lfir) for each 
source. We have used two different methods to estimate FIR lumi- 
nosities based on: converting the Spitzer 24 flux to (Lfir) and, 
interpolating from the measured 850 fim flux density. Although 
Lfir estimates of high-z objects carry uncertainties, depending 
on the method chosen, it is reassuring that the methods presently 
considered, yield within errors, similar results. In what follows we 
discuss each method in detail. 

At the redshifts of our type-2 QSOs the Spitzer 24/im band 
corresponds to rest-frame mid-IR (6-8 /xm depending on redshift) 
which is, in general, used to infer the total thermal (8-1000 /itm) 
Lfir (e.g. Chary & Elbaz 2002). A number of authors (e.g. 
Spinoglio et al. 1995, Rush et al. 1993) have shown that, especially 
for AGN, the MIR (12/im) luminosity is a good proxy for the FIR 
luminosity. Following AH06 we estimate the 12/im-to-FIR lumi- 
nosity conversion factor according to the two templates (Arp220 
and IR22491) we used to fit the SEDs of our type-2 QSOs. We 
derive the 12 fim rest-frame luminosity by extrapolating from the 
observed 24 /im flux. The resulting Lfir values are listed in Ta- 
ble 3. All our type-2 QSOs are highly luminous, with luminosities 
formally in the "ultra-luminous" class (Lfir > lO^^L©). 

In section 3 we showed that the far-infrared/ submm part of the 
SED of type-2 QSOs can be well represented by the two ULIRG 
templates Arp220 and IR2249 1 . In both of these template galaxies 
submm emission is emitted by dust heated by young stars (rather 
than the AGN continuum). Under the assvunption that the submm 
emission in type-2 QSOs is also due to star-forming activity we 
can use the 850 /xm fluxes to estimate Lfir (by scaling them with 
those of the templates used to fit each of the type-2 QSOs). For 
IR22491 we use log LFiR=12.n and for Arp220 log Lfir=12AS 



(e.g. Rigopoulou et al. 1999). Lfir values estimated in this way 
are of course subjected to uncertainties depending on the template 
of choice. We note however, that the choice of template only af- 
fects the normalization of the derived luminosities; its variation 
as a function of redshift over the range covered by our sources 
(1.0<z<3.5) is small for both templates considered, because the 
strong negative K-correction at 850 fira cancels out the effect of 
cosmological dimming (e.g. Blain & Longair 1993). 

Additionally, for those sources with radio detections we have 
used the well-established radio-far-infrarcd correlation to estimate 
Lfir with values in good agreement with the previously discussed 
methods. We list the computed Lfir values in Table 3. 

Using the Lfir values in Table 3 and following the prescrip- 
tion by Kennicutt (1998) we estimate SFR rates by 

SFR = LF/i{/(5.8 X IO^Lq) {MQ/yr). 

Using the observed 850 jim fluxes the dust mass is then given 

by: 

Mdust = (1 + Z) ^S850D'l/krestB(Urest,Ta) 

where z is the redshift, D_l is the luminosity distance, krest is 
the rest frequency absorption coefficient and B(Vr-est,Td is the rest 
frequency value of the Planck function from dust grains radiating 
at temperature T^. Using krest=0.l5 (e.g. Scott et al. 2002) and as- 
suming optically thin thermal emission the derived Mdust are listed 
in Table 3. 

Finally, following Marconi et al. (2004) we estimate bolomet- 
ric luminosities based on rest-frame absorption corrected Lx- As- 
suming that Lboi = LEdd, where LEdd is the Eddington luminosity 
we report (Table 3) lower limits on the black hole masses. 



6 ORIGIN OF THE SUBMILLIMETRE EMISSION 

A different approach to the relative importance of AGN and star 
formation activity is provided by comparing the submillimetre (850 
/xm) to X-ray (2 keV) spectral slope (asx) of our type-2 QSOs to 
those of template nearby galaxies. We have, therefore, calculated 
the asx values for our sample galaxies using: 
asx = -log[i^^ - 0.18] x 0.163 
where f2kev is the observed flux density at 2 keV (keV cm~^ 
keV" '^), calculated from the full-band flux using the estimated pho- 
ton index, and fgso is the observed flux density at 850 /im (mJy). 
Figure 4 shows the estimated asx values as a function of redshift 
for our sample t3T)e-2 QSOs. In the same plot we also show tem- 
plate SEDs for: an unabsorbed quasar, quasars with varying degrees 
of absorption (Nh = 10^^, 10^" cm"^ and Compton-tiiick with 1% 
of the nuclear emission scattered), a mean starburst and an Arp220- 
like. 

The templates used here are discussed in more details in Al- 
maini et al. (2003). The type-2 QSOs yield values ranging be- 
tween 1.41 (for CDFS-901) to 1.17 (for CDFS-45). Larger values 
of Oisx indicate stronger submillimetre emission relative to the X- 
ray emission. The asx indices are clearly incompatible with un- 
absorbed AGN. The asx index for CDFS901, CDFS54, CDFS51 
and CDFS112 are compatible with the mean starburst and Arp220 
templates. It is of interest to note that of these four sources three 
are best fit by the Arp220 template (see discussion in section 3). 
The remaining sources CDFS263, CDFS 27 and CDFS45, have in- 
dices that are close to the spectral index of a Compton-thick AGN 
in which only 1% of the nuclear emission is seen through scatter- 
ing ( the properties of CDFS263 have already been discussed in 
Mainieri et al. 2005). For comparison we report in Figure 4 the 
spectral indices calculated for the submillimetre sources studied by 
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Figure 3. Colour-colour plots involving IRAC (left) and IRAC and MIPS bands (right, adapted from Lacy et al. 2004). Symbols are as follow: filled circles: 
type-2 QSOs, filled squares: warm SCUBA sources, open squares: cold SCUBA sources (mid-infrared Spitzer data for SCUBA sources from Egami et al. 
2004 and Ashby et al. 2006). Those type 2 QSOs that follow the Arp220 template SED (see discussion in section 3) are marked with a squared circle. The 
wedged-shaped (plot 3a) and the hatched areas (plot 3b) denote the regions occupied by AGN. 



Table 3. FIR Luminosities, SFR, Dust and Black Hole Masses 



ID 


logLi^/j{(12/im)^ 
L© 


logLi^/j{(850/im) 
L© 


SFR2 
MQ/yr 


xlO^M© 


X 10^ M© 


54 


12.58 


12.78 


670-1165 


4.9 


0.5 


45 


12.49 


12.72:^ 


540-905 


4.1 


1.4 


263 


12.55 


12.64 


620-1190 


4.2 


1.9 


27 


12.34 


12.62: 


380-720 


5.1 


0.7 


112 


12.67 


12.88 


710-1370 


5.5 


0.8 


901 


12.60 


12.87 


690-1460 


4.8 


0.6 


51 


12.40: 


12.79: 


440-920 


3.9 


1.6 



: luminosities based on 24 fj,m fluxes, see text 
^: range of SFR based on different Lpm 
^: (:) denotes upper limit 
black hole masses are lower limits 



Alexander et al. (2005). These sources show a wide range in in- 
dices similar to what we reported for our type 2 QSOs. In terms of 
absolute asx value, the four type-2 QSOs we discussed before are 
indistinguishable from the submillimetre sources. 

Another approach to assess the relative importance of the two 
components (AGN and starburst) is to examine the x-ray-to-far- 
infrared luminosity ratio. While this ratio is in principle quite sim- 
ilar to the asx spectral index, it is independent of the X-ray spec- 
tral slope (which one has to assume when converting broad band 
to monochromatic flux densities). In Figure 5 we plot the Lx vs 
Lfir ratio for the present sample of type-2 QSOS (with values 
taken from Tables 2 and 3) and Uterature galaxies which also host 
obscured AGN. Such a comparison is instructive with the caveat 
that the Lf/« estimates may originate from a variety of measure- 
ments (e.g mid-infrared, submm, etc). In the same plot we show a 
mean "starburst" L^j Lfir ratio (adapted from the work of David, 
Jones and Forman 1992) and the same ratio for quasars (from Elvis 
et al. 1994). It is difficult to define the same mean Lx/Lfir ra- 
tio for AGN since they show such a huge dispersion in their values. 
Additionally, we have marked the location of Arp220, Mrk23 1 and 
NGC 6240. The submillimetre sources (in particular those whose 
X-rays originate from an AGN, Alexander et al. 2005) display val- 
ues very similar to those of the present type-2 QSOs. For further 



comparison we have also marked the location of the SDSS opti- 
cally selected type-2 QSOs (Zakamska et al. 2008). These show 
a wide spread in their L^/ Lfir, spanning from values as low as 
those displayed by starbursts all the way to the quasar values. None 
of the type-2 QSOs have L^JLfir values larger than those of the 
quasars presented in Elvis et al. (1994). This finding implies that 
type-2 QSOs are likely to experience some level of star-formation 
activity which could be responsible for the FIR and submillimetre 
emission. In turn, such star-forming processes could well contribute 
to the total bolometric output although, at this stage we cannot ac- 
curately attribute a specific fraction. Within our very small sample 
we do, however, see a range in the strength of the starburst which 
may reflect the evolutionary stage of each object. In the next sec- 
tion we investigate in more detail a possible evolutionary scheme 
for these sources. 



7 THE EVOLUTION OF TYPE-2 QSOS 

It has been suggested that the early phases of QSO evolution are 
characterised by substantial absorption (e.g. Fabian 1999). In such 
models the main obscured growth phase of the QSO coincides with 
the formation of the galaxy spheroid. Once this phase is completed 
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QSO Njj=102* 




•St 



redshift 

Figure 4. SubmiUimetre to X-ray spectral index {asx as a function of redshift. The curves show the expected asx values for a set of SEDs (QSOs with 
various degrees of obscuration, mean Starburst, Arp220). The values for the present type 2 QSOs are denoted by solid squares. Filled circles correspond to 
SCUBA galaxies from Alexander et al. (2005). 




1040 10*' 10*2 10*3 10** 10*5 1046 1047 

LJ2.0- 10 keV](erg/s) 

Figure 5. Rest-frame far-infrared vs. absorption-corrected [2-10] keV luminosity for our type-2 QSOs (filled squares) and other obscured AGN from the 
literature. SDSS optically selected type-2 QSOs (open squares) from Zakamska et al. (2008), SubmilUmetre galaxies (filled circles) from Alexander et al. 
(2005), local Universe "template" objects, Arp 220, Mrk 23 1 and NGC 6240 (filled triangles). The diagonal lines show ratios of X-ray to far-infrared luminosity 
for Starbursts (adapted from the work of David, Jones and Forman 1992, the X-ray data have been converted to the [2-10] band used here assuming F = 2.0) 
and quasars (adapted from Elvis et al. 1994). 



then the QSO blows away the obscuring material and starts its lu- 
minous unobscured phase (e.g. Silk & Rees 1998). Such models 
deviate from the traditional view of the "unified scheme" for local 
AGN (e.g. Antonucci et al. 1993) whereby the absorbing material 
is present in all objects in the form of a toroidal structure. Based on 
simulations of galaxy mergers that include black holes (e.g. di IVlat- 
teo et al. 2005; Springel et al. 2005), Hopkins et al. (2005) deter- 
mine the intrinsic lifetime of the strong accretion phase and, taking 
into account obscuration, infer a lifetime for the QSO. They find 



that the optical QSO phase is rather short compared to the intrinsic 
lifetime of the quasar. Additionally, their model suggests the pres- 
ence of an obscured phase (type 2 QSOs) that occurs irmnediately 
before the dust is completely blown away. 

By contrasting the X-ray/ submillimetre properties of an X- 
ray absorbed, Compton-thin QSO sample with that of a matched 
unabsorbed sample of QSOs, Page et al. (2004) and Stevens et 
al. (2005) provided observational evidence for the last two stages 
of the evolutionary scenario described above. The X-ray absorbed 
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QSOs are observed during the transition phase, where the QSO is 
still actively forming stars (formation of the galactic bulge) while 
at the same time accretes material that assists the formation of the 
central black hole. Once the black hole has reached a mass > few 
X 1O*M0 the QSO enters its unobsucred luminous phase, slowly 
switching off the star-forming process and eventually evolving into 
relaxed elliptical galaxies (e.g. Kukula et al. 2001, Dunlop et al. 
2003). 

With bolometric lurrunosities in excess of Ltoi > 
10*® erg's"^, obscuring colvmrn densities Nif >10^^cm~^, the 
present type-2 QSO sample can be used to further investigate the 
QSO evolutionary scheme. In section 3 we compared the SEDs of 
type-2 QSOs with those of local ULIRG templates 1R22491 and 
Arp220. Local ULIRGs are characterised by extremely high bolo- 
metric luminosities and similarly large colunm densities. If we con- 
vert the estimated Av of ULIRGs (e.g. from Genzel et al. 1998) 
to expected Nj/ that results in N_ff values in excess of lO^'^cm"^. 
Given the high obscuration, the present (albeit limited at the mo- 
ment) detections in the submillimetre, large dust masses (and con- 
sequently gas) and high redshifts, the present sample of type-2 
QSOs could be candidates for the initial highly obscured phase of 
QSO evolution. A pre-requisite for any QSO to be "caught" in this 
phase is of course submillimetre emission. During this phase the 
black hole must be accreting rapidly in order to achieve a mass 
of ~1O*M0 and thus should also reach X-ray Itmiinosities in ex- 
cess of 10'*'' ergs^'. And likewise the local ULIRGs, the proper- 
ties of QSOs in this initial phase of evolution are not homoge- 
neous in terms of Ituninosity, submillimetre brightness, even the 
starburst/AGN contribution. 

Is there a link between type-2 QSOs and submm-bright galax- 
ies? In section 4 we compared the mid-infrared properties of the 
two samples and found that they show certain similarities. Addi- 
tional evidence was provided in Figures 4 and 5 through the com- 
parison of the asx indices and the x-ray-to-far-iirfrared luminosity. 
For more than half of the type-2 QSOs their asx index is identical 
to that measured for submillimetre-bright objects from the Alexan- 
der et al. (2005) sample. It is thus quite possible that a fraction of 
the submm-luminous sources are in fact caught during this initial 
highly absorbed and actively star-forming phase in the QSO evo- 
lution. The evidence provided by Alexander et al. (2005) on the 
existence of black holes in a fraction of the submillimetre-bright 
population gives further support to this claim. Likewise, we have 
found that the Lx/Lf/a ratio of the present type-2 QSOs is smaller 
than that found for the well studied quasars of Elvis et al. (or objects 
that are purely AGN dominated). 

We have presented some evidence for the initial stages of 
QSO evolution by reporting measurements of submillimetre emis- 
sion from a sample of highly absorbed type-2 QSOs discovered 
through deep X-ray surveys. Two out of the eight targets have firm 
detections {S/N >4) and a further three objects are marginally de- 
tected. These submillimetre detections coupled with their overall 
similarities to local ULIRGs suggest that these QSOs may 

be experiencing co-eval growth of their spheroid component while 
at the same time maintaining an accretion rate high enough to form 
a black hole of Mbh ~1O*M0. Once this is achieved, the QSO 
will blow away the obscuring material revealing its luminous QSO 
and following a passive evolution into the present-day ellipticals. 
Future sensitive submillimetre observations are clearly needed to 
further strengthen this claim. Deep sensitive surveys at 250, 320 
and 520 /um with HERSCHEL and SCUBA -2 (450 and 850 /ttm) 
will be crucial in delineating the shape of the far-infrared/submm 



part of the SED of type-2 QSOs and provide firm evidence of their 
submillimetre emission. 
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